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Abstract 

In this paper the effect of the post-Newtonian gravitomagnetic force 
on the mean longitudes I of a pair of counter-rotating Earth artificial 
satellites following almost identical circular equatorial orbits is investi- 
gated and the possibility of measuring it is examined. The observable 
is the difference of the times required for I to pass from to 2tt for both 
senses of motion. Such a gravitomagnetic time shift, which is indepen- 
dent of the orbital parameters of the satellites, amounts to 5x10^^ s 
for the Earth; it is cumulative and should be measurable after a suf- 
ficiently high number of revolutions. The major limiting factors are 
the unavoidable imperfect cancellation of the Keplerian periods, which 
yields a constraint of 10~^ cm in knowing the difference between the 
semimajor axes a of the satellites, and the difference / of the inclina- 
tions i of the orbital planes which, for i ^ 0.01°, should be less than 
0.006°. A pair of spacecraft endowed with a sophisticated intersatellite 
tracking apparatus and drag- free control down to 10^^ cm s^^ Hz^s 
level might allow to meet the stringent requirements posed by such a 
mission. 
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Figure 1: Orbital geometry for a motion around a central mass. Here L 
denotes the orbital angular momentum of the particle of mass m, J is the 
proper angular momentum of the central mass M, 11 denotes the pericentre 
position, / is the true anomaly of m, which is counted from 11, Vt, uj and i 
are the longitude of the ascending node, the argument of pericentre and the 
inclination of the orbit with respect to the inertial frame {x,y,z} and the 
azimuthal angle (f) is the right ascension counted from the x axis. 

1 Introduction 

The orbital path of a test particle freely falling in the gravitational field of 
a central body of mass M and angular momentum J is affected by various 
post-Newtonian general relativistic effects of order 0(c~^)UI21 like the Ein- 
stein gravitoelectric secular precession of the argument of pericentre^ uj jS] 
and the Lense-Thirring gravitomagnetic secular precessions of the longitude 
of the ascending node and the argument of the pericentre uj :.6,_ . The Ein- 
stein precession is now known at a 10~'* level of relative accuracy, e.g., from 
Solar System planetary motions by analyzing more than 280,000 position 
observations (1913-2003) of different types including radiometric observa- 
tions of planets and spacecrafts, CCD astrometric observations of the outer 
planets and their satellites, meridian transits and photographic observations 
[7] . On the contrary, to date there are not yet clear and undisputable direct 
tests of the various effects induced by the gravitomagnetic field. In April 
2004 the Gravity Probe B (GP-B) mission has been launched. Its goal 

^For an explanation of the Keplerian orbital elements of the orbit of a test particle see, 
e.g., |1| and Figure 
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is to measure, among other things, the gravitomagnetic precessions of the 
spins of four superconducting gyroscopes [HI carried onboard with a claimed 
accuracy of 1% or better. The Lense-Thirring effect on the orbits of the 
LAGEOS and LAGEOS II Earth's artificial satellites has been experimen- 
tally checked for the first time some years ago with a claimed accuracy of 
the order of 20% TTT, but, at present, there is not yet a full consensus on 
the real accuracy obtained in that test jllj . 

Recently, it has been shown that also the orbital period of a test par- 
ticle is affected by the post-Newtonian gravitoelectromagnetic forces. The 
gravitomagnetic correction to the period of the (equatorial) right ascension^ 
a, for e = 0, i = 0°, where e and i are the eccentricity and the inclination, 
respectively, of the test particle's orbit, has been worked out in ^^l- The 
case of more general orbits has been treated in ^Ij (e = 0, i 7^ 0°) and ^Hl 
(e 7^ 0, i 7^ 0°). The gravitoelectric correction has been treated in _T3]. Let 
us consider a pair of counter-orbiting satellites, conventionally denoted as 
(+) and (— ), following identical circular and equatorial orbits along opposite 
directions around a central spinning body of mass M and proper angular 
momentum J; while the gravitoelectric correction to the Keplerian period 
has the same sign for both the directions of motion, similar to the classi- 
cal perturbing corrections of gravitational origin (oblateness of the central 
mass, tides, N-body interactions), the gravitomagnetic one changes its sign 
if the motion is reversed. This allows, at least in principle, to single out the 
latter one by measuring the difference of the periods of the counter-orbiting 
satellites. Some preliminary error analyses can be found in |16j . 

As we will show here, it turns out that also the mean longitudes^ I = 
M. + 0, + uj are affected, among other perturbations, by the post-Newtonian 
general relativistic gravitoelectromagnetic forces. Also in this case the grav- 
itomagnetic correction is sensitive to the direction of motion, contrary to 
the gravitoelectric one. 

In this paper we will investigate the possibility of singling out the gravit- 
omagnetic effect on the mean longitudes by means of a suitable space-based 
mission in the Earth space environment. We will mainly deal with a number 
of competing classical effects which could represent a source of systematic 

^It is nothing but the azimuth angle </!> of a spherical coordinate system in a frame 
whose origin is in the center of mass of Earth, the {x, y} plane coincides with the Earth 
equatorial plane and the x axis points toward the Vernal Equinox Y (see Figure 0. In 
satellite dynamics it is one of the direct observable quantities |12| . 

"^M is the mean anomaly of the test particle's orbit. I, along with the Keplerian orbital 
elements, comes out from the machinery of the data reduction process performed by the 
orbit determination softwares like GEODYN II and UTOPIA. 
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errors and bias. No measurements error budget is carried out. 



2 The Gravitomagnetic Effect on the Mean Lon- 
gitude 

From the Gauss perturbative equations it turns out that the rate equa- 
tion for /, for small but finite values of e and i, is 

~ n Ar -+(-)— + 



dt na \aJ \2 J dt \2 J dt 

where a is the satellite's semimajor axis and n = {GMa~^)2 is the Keplerian 
mean motion; Ar is the radial component of the perturbing acceleration 
apert- The currently available technologies allows to insert Earth artificial 
satellites in orbits with^ e < 10^^ and i ~ 0.01°; then, the third and fourth 
terms in eq.(^, which are proportional to i and e because dQ/dt oc 1/ sini 
and duj/dt oc 1/e, can be neglected and the equation 

r^n Ar (2) 

dt na 

can be used instead of eq.(p3) to a good level of approximation. 

According to ^H], the radial components of the post-Newtonian general 
relativistic gravitoelectromagnetic accelerations, for generic orbits around a 
central spinning body, are 

^(CE) ^ (GMni + e cos /)^ ^3 ^ ^ _ ^ 

^(GM) ^ ^"^'^^°^\ (l + ecos/)^ (4) 
c2a2(l — e^) 2 

where / is the true anomaly. They are induced by the post-Newtonian 



This value holds for initial orbital injection. Later, various perturbations, mainly due 
to the terrestrial gravitational field, afi^ect e which can reach values up to 5 x 10""^. 
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general relativistic gravitoelectromagnetic fields^ 

^ = -^(^-^ Jr--^(r.v)v, (6) 

b(gm) = _^[j_3(j.f)f]. (7) 

Prom eq. ((nj-eq. it can be noted that, while the gravitoelectric acceler- 
ation is insensitive to the sense of motion of the test particle along its orbit, 
it is not so for the gravitomagnetic acceleration due to its dependence on 
cosi. Indeed, according to ^H] 

■ • (^'+^')^ 
smz = phj ' 

cos* = TTly (9) 
I h| 

where h is the orbital angular momentum per unit mass whose components 
change sign when v —v. By reversing the sign of the velocity vector one 
obtains (+ and — denote the pro- and retrograde orbits, respectively) 

sin = sini^~'"\ (10) 

cosi*-^-* = — cosi^'*'^ (11) 

from which it follows 

i(-) = 180° -iW. (12) 

By inserting eq.Q-eq.lji} in eq.©, it can be obtained for the time Pi 
required to / for passing from to 2it 

p(±) p(0) + p{GE) ^ (GM) ^2n^ U^jGMro)'^ ^ 
Ea.((T^ yields for identical orbits followed in opposite directions. 



AF, := Pr - Pr = (14) 

which amounts to 5 x 10~^ s for Earth and is four times larger than the 
corresponding effect for the right ascension. 

^In the weak-field and slow-motion approximation of the General Theory of Relativity 
the equations of motion of a test particle are 

J = -Eg-E('3^)-2:^ xB(G^), (5) 
where Eg = GMr/v^ is the usual Newtonian monopole term. 
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3 Sources of Systematic Errors 



Eg . (1141) would be valid only if the Keplerian periods (and the various classical 
and post-Newtonian gravitoelectric perturbative corrections to them) of the 
two satellites were exactly equal; this condition, however, cannot be achieved 
due to the unavoidable orbital injection errors. Then, ea. (ll4() has also to 
account for the difference induced in the Keplerian periods and the various 
classical and post-Newtonian gravitoelectric corrections to them, e.g., by 
the difference d in the semimajor axes of the two satellites. The differences 
in i, which do not affect the Keplerian periods and the post-Newtonian 
gravitoelectric correction, would have, instead, an impact on the classical 
perturbative terms which cannot be neglected. 

In general, the difference between the mean longitude periods of the two 
counter-rotating satellites can be written as 

APi = AP(°) + APf 1^"" + AP/™^ + AP/^^^ . (15) 

Over many orbital revolutions, say N, the accuracy in determining the grav- 
itomagnetic time shift becomes 

^[Ap(GM)] < '^(AP)""P ^ ^[^p(o)] ^ ^[APfi'^^^] + 6[APp^], (16) 

where 5(APi)'^^P is the experimental error in the difference of the obtained 
P^^^ over revolutions. After a sufficiently high number of orbital revo- 
lutions it should be possible to make the term d^APiY^"^ /N smaller than 
^p^(GM) order to get an estimate, let us calculate the angular shift 
corresponding to AP^'''^'^'' over an orbital revolution by using 

For^ ro = 25498 km and = 4.05200895378 x 10"^ s, I = 27r the grav- 
itomagnetic shift amounts to A^ = 1 x 10^^ milliarcseconds (mas in the 
following; 1 mas = 4.8xl0~^ rad), over an angular span. Since the accu- 
racy with which the orientation of the axes of the International Terrestrial 

®In order just to fix the ideas, we shall consider the semimajor axis of the existing SLR 
ETALON satellites. We will see that high altitudes contribute to reduce the impact of 
many aliasing perturbations. 
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Reference System (ITRS) is of the order^ of 3 mas, after 300 orbital revolu- 
tions or 140 days, the gravitomagnetic time shift would reach this sensitivity 
cutoff. 



3.1 The Impact of the Imperfect Cancellation of the Keple- 
rian Periods 

As we will see, the major limiting factor in measuring the gravitomagnetic 
time shift of interest is the difference of the Keplerian orbital periods 



Ap(o) 



(ro + d)^ 



{GM)2 



{GMy. 



(18) 



where tq represents the nominal value of the semimajor axis of the two 
satellites. The shift of ea. (|18|) cannot be made smaller than ea. H14() . in 
practice, by choosing suitably the orbital geometry of the satellites. Indeed, 
from ea. (|18l) follows that the relation dr^^ < 1.059 x 10^ cma should be 
fulfilled; for ro ~ 10^ cm it would imply d ~ 3x 10 cm. Then, in order to be 
able to measure the relativistic effect of interest, which accumulates during 
the orbital revolution, the difference should be subtracted from the data 
provided that its error 5[APW], which is present at every orbital revolution 
and is due to the uncertainties in the Earth GM, ro and d, is smaller than 
the gravitomagnetic time shift. This error is given by 



SirdrQ 



2(GM)2 



6iGM) + 



Snd 



2{GMro) 



X {6ro) + 



37rrQ 



{5d) 



and by assuming ro = 25498 km, {droY^'^ 



(GM)- 

(19) 

1 cm and 6{GM) = 8 x 10^^ 



cm" s 



EOl, eq.dHI) yields 



(5[AP(°)] < (2.8 X 10"^^ cm~^ s) x d+{2.3 x 10" 



cm 



s) X {6d). (20) 



Eg.ipn]) tells us that the error due to the uncertainty in^ GM and, espe- 
cially, ro is less than the gravitomagnetic effect, while 5d should be at the 

^Its realization, the International Terrestrial Reference Frame (ITRF), is obtained by 
estimates of the coordinates and velocities of a set of observing stations on the Earth |2U| . 
*See on the Internet http://hpiers.obspm.fr/iers/bul/bulb/explanatory.html 
^Note that the uncertainty in the Earth's GM is a not a totally negligible limiting 
factor. Indeed, for an intersatellite separation of a few kilometers the induced bias is 
of the order of 10~* s. If only one satellite was considered the systematic error in the 
Keplerian period due to GM ^(a/GM)^ 5{GM) would amount to 4 x 10^^ s, as pointed 
out in the preliminary evaluations of |22|. 
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level of 2 X 10^^ cm. It seems to be impossible to meet this very stringent 
requirement with the current SLR technology due to many measurement 
errors (station errors, random errors in precession, nutation and Earth rota- 
tion, observation errors). However, an intersatellite tracking approach could 
yield better results. A level of of the order of 10~^ cm is currently 

available with the K-band Ranging (KBR) intersatellite tracking technology 
used for the GRACE mission |2J|. On the other hand, it must be noted that 
it would not be easy to use this approach for satellites following opposite 
directions. The GRACE twins revolve in the same sense. 

In regard to ( ,5 systematic^ thinking about a pair of completely passive, 
spherical LAGEOS-like satellites, it should be pointed out that for orbits 
with e ~ i ~ 0, for which ea. ()14|) holds, many non-gravitational perturba- 
tions^^ affecting the semimajor axis vanish |22j : the remaining ones could be 
strongly constrained by constructing the two satellites very carefully with 
regard to their geometrical and physical properties. Some figures will be 
helpful. We will adopt the physical properties of the existing LAGEOS 
satellites. The atmospheric drag induces a decrease in the semimajor axis 
of Aa = —a? qCdS/iti over one orbital revolution, where Cd is the satellite 
drag coefficient, S/m is the satellite area-to-mass ratio and q is the atmo- 
spheric density at the satellite altitude. For Co ~ 4.9, S/m = 7 x 10"'^ cm^ 
g~^ and Q = 8.4 x 10~^^ g cm~^ (estimate of the atmospheric neutral den- 
sity at LAGEOS altitude; for tq = 25498 km it should be smaller, of course) 
the decrease in the semimajor axis amounts to 4xl0~^ cm over one orbital 
revolution. The impact of direct solar radiation pressure on the semimajor 
axis of a circular orbit vanishes, also when eclipses effects are accounted for. 
The nominal Poynting-Robertson effect on the semimajor axis of a circu- 
lar, equatorial orbit (ro = 25498 km) is of the order of 10~^ cm over one 
orbital revolution. However, it is linearly proportional to the satellite re- 
flectivity coefficient Cn for which a 0.5-1% mismodelling can be assumed. 
The semimajor axis of a circular orbit is neither affected neither the albedo 
over one orbital revolution. The same holds also for the IR Earth radiation 
pressure and the solar Yarkovsky-Schach effect (by neglecting the eclipses 
effects). The Earth Yarkovsky-Rubincam effect would affect the semimajor 
axis of an orbit with i ~ 0° by means of a nominal secular trend less than 
1 X 10~^ cm (ro = 25498 km) over one orbital revolution^^ . However, a 20- 
25% mismodelling can be assumed on it. The direct effect of the terrestrial 

^"^On the semimajor axis there are no long-term gravitational perturbations. 
^^It would vanish if the satellite spin axis was aligned with the z axis of an Earth- 
centered inertial frame. 
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magnetic field on the semimajor axis of a circular orbit vanishes over one 
orbital revolution. 

Another possibility could be the use of more complex (and more expen- 
sive) spacecraft endowed with a drag-free apparatus; it would be helpful in 
suitably constraining ((^7.Q)systematic ^^^^^ especially, ((^c?) systematic^ Indeed, for 

orbits with semimajor axis of the order of 10^ cm the unperturbed Keplerian 
period is of the order of P^^^ ~ 4 x 10^ s. If the maximum admissible error 
in knowing the position of the satellites is of the order of 2 x 10~^ cm then 
the maximum disturbing acceleration that would not mask the gravitomag- 
netic clock effect, over one orbital revolution, is 6d/[P^^^'^ ~ 1 x 10^^^ cm 
s~^. But the drag-free technologies currently under development for LISA 
[23) and OPTIS missions should allow to cancel out the accelerations 
of non-gravitational origin down to 3 x 10^^^ cm s~^ Hz~2 and 10~^^ cm 
s^^ Hz^2, respectively. In our case, for tq = 25498 km and an orbital fre- 
quency of 2.4 X 10~^ Hz, a 1 X 10~^ cm s~^ Hz" 2 drag-free level would be 
sufficient. More precisely, from the Gauss perturbative rate equation of the 
semimajor axis it turns out that, for circular orbits, the effects induced on 
''0 by periodically time-varying non-gravitational along-track accelerations 
with frequencies multiple of the orbital frequency vanish over one orbital 
revolution. For constant along-track accelerations At the situation is differ- 
ent: over one orbital revolution (corresponding to a frequency of the order 
of 10^^ Hz) their effect on tq would amount to 

In order to make Sr^'^^ = 2 x 10~^ cm the maximum value for At would 
be 3.8 x 10~^^ cm s~^ in a frequency range with 10"^ Hz as upper limit. 
It is a difficult but not impossible limit to be obtained with the drag-free 
technologies which are currently under development. 



3.2 The Impact of the Imperfect Cancellation of the Post- 
Newtonian Gravitoelectric Periods 

The perturbative correction to Pi induced by the post-Newtonian gravito- 
electric acceleration, given in ea. H14() . amounts to 

^(GE) ^ 12^(G'Mro)i ^ 
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For an Earth high-orbit satehite it is of the order of 10 ^ s. The orbital 
injection errors in tq would yield 

^ ^(GE) evr ( GM\^ , , , 

which amounts to 10~^ s for tq = 25498 km and d = 5 km. This result 
justifies, a posteriori, our choice of neglecting the small corrections to ea .(12211 
induced by the small, but finite, eccentricity of the orbit. 

3.3 The Impact of the Classical Gravitational Perturbations 

In this section the perturbations of gravitational origin on the mean longi- 
tude / and their impact on the measurement of l^P^'^^'' are investigated. 

3.3.1 The Geopotential 

Since we are interested in effects which are averaged over many orbital rev- 
olutions, only the zonal harmonics of geopotential, which induces secular 
perturbations on Vt, to and Ai, will be considered. In the following the 
Gaussian approach together with eq.© will be used. 
Let us write the zonal part of geopotential as 

$zon ._ ^(0) ^ ^^zon ^ ^ _ (^^f ^ '^'^(^08 0). (24) 

^ e=2 ^ 

According to (6.98b) of the radial component of the perturbing accel- 
eration, in spherical coordinates, for an even zonal harmonic of degree i, is 

<^"^^^T^M^ + l)^^<A(cos., ,25, 

It is easy to see that the odd zonal harmonics of geopotential do not affect 
the mean longitude of a satellite in an equatorial orbit. Indeed, it is well 
known that the odd-degree Legendre polynomials vanish for 9 = 90° j2()j . 
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The perturbing radial accelerations due to the even zonal harmonics are 



4'^^ = (26) 



^ (27) 

' (28) 

_ 315GMfi|j8 

Inserting ea.(PB |) -ea. H29|) in eq.@ yields 

„bi _ QttR%J2 157ri?|j4 357rii^J6 SlSvri^lJs 



lhGMR%Ji 



{GMro)2 (GMr5)2 4{GMr^)2 'i2{GMrf 



= -8.2341063 s - 1.9266 x lO^^ s - 2.34 x 10^^ s + 5 x lO""^ s + ... (30) 

for an orbit radius ro = 25498 km. This shows that, for such a radius, just 
the first four even zonal harmonics are of the same order of magnitude of 
P/^^) ~ 10-7 s. 

Let us now evaluate the impact of the even zonal harmonics the geopoten- 
tial on the measurement of /S.P^'^^^ for a pair of counter-orbiting satellites 
whose orbits differ by a small amount d <C in radius. With 

, - 7 ~ , , n + 2 > I3ij 

(r-o)2 {ro + d)2 2{ro)~ 

the difference in the perturbing terms due to the even zonal harmonics are, 
for d = 5 km and tq = 25498 km 

AP/^^) = = _8.068 X 10-S, (32) 

^p(Ji) ^ ^4 ^ _g ^ ;^Q-7 33 

2{GMrl)2 

( Ja) 3157ri?!l Jfifi o 
^p(-J6) ^ = -2 X 10"^ s, (34) 







8(GMrii 

40957ri?|J8C^ 
64(GMr^5)| 



Ap/^«) = ®_ ? = 6 X 10-1° s. (35) 
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This means that only AP^^^ and AP^^'^ are relevant. The uncertainty in 
the knowledge of the Earth GM and of the even zonal harmonics yield 

6[ApI^''^]gm = ^^^^^4 X 5(G'M) ~ 8 X 10-13 s, (36) 
2(GMro)t 

5[APi'%, = -^^^ X ~ 3 X 10-" s, (37) 

6[APy'^]GM = x6{GM)r.9xW-^^s, (38) 

4(GM)2(ro)2 

S[AP['%, = ^^^^ X (5J4) ~ 3 X 10-11 s, (39) 



For 6{GM) and 5J£ the values of the lERS convention jJOj and of the 
EIGEN-2 Earth gravity model j27j . respectively, have been used. 

Another source of error is represented by the uncertainty in the knowl- 
edge of the separation d between the two orbits and of their radius rg . Their 
impact is, by assuming 6d ~ 6rQ ~ 1 cm 

= ^"^e^^^ x(M -4x10-13 s, (40) 
2(GM)2(ro)- 



2 



5[AP/^^)], = ^^x{6d)^lxW''s, (41) 
(GMrj])" 



I 2 



= X (5ro) ~ 1 X 10-1^ s, (42) 

4(GM)2(ro)- 



2 



5[AP/^*)], = X (M) ~ 2 X 10-12 s, (43) 

2(GMr7)2 



These results clearly show that, for the same semimajor axis of, say, 
the ETALON SLR satellites and for not too stringent requirements on the 
orbital injection errors for the separation between the two orbits, the impact 
of the geopotential on the measurement of AP^'^^^ is negligible. 

Finally, let us consider what could be the impact of the secular vari- 
ations of the even zonal harmonics J^. For the first even zonal harmonic 
the effective jf ~ ^2 0.371^4 + 0.079^6 + O.OOGJs - 0.003jio..., whose 
magnitude is of the order of (—2.6 =b 0.3) x lO-n yr-i, can be considered. 

/ rcff \ 

The insertion of such value in eq. dSH) yields APf ^ ^ = i x 10"" s over one 
year, so that it can be concluded that its effect can be safely neglected. 
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3.3.2 The Tides 



Another source of potential systematic error is represented by the solid Earth 
and ocean tides j29j . 

For a constituent of given frequency f of the solid Earth tidal spectrum, 
which is the more effective in perturbing the Earth satellites'orbits, the 
perturbation of degree i and order m induced on I is 

^ +00 

j^=n+A,^{GM)'2R^-\-^'-^h{l+l)Y, FimpGipgk^S.Hr cos ^umpg, 

p=0 q=—oo 

(44) 

where Aim = [(2^ + - m)l/{iir){£ + m)!]2, Fimpii) and G£pg{e) are 
the inclination and eccentricity functions, respectively 0, are the Love 
numbers, are the tidal heights and jumpg is built up with the satellite's 
orbital elements fi, uj and M, the lunisolar variables and the lag angle 6fim 
pn] . Note that the dependence on Fimp and Gipg is the same also for the 
ocean tidal perturbations. 

Prom an inspection of the explicit expressions of the inclination and 
eccentricity functions jlj and from the condition i — 2p + q = 0, which must 
be fulfilled for the perturbations averaged over an orbital revolution, it turns 
out that the i = 3 part of the entire tidal spectrum does not affect the mean 
longitude of a satellite with i = e = 0. With regard to the i = 2 constituents, 
only the long-period zonal {m = 0) tides induce non vanishing perturbations 
on Pi for circular and equatorial orbits. Among them, the most powerful is 
by far the 18.6-year lunar tide. Its effect on Pi is given by 

p(18.6-yr) _ 3(5vr)li?e 4?(f)^2°(f) ^ 

M — 71 cos 7f 2010- \^^) 

[GMro)2 

For vq = 25498 km it amounts to 3x10^^ s. The difference in the orbits of 
the two satellites would yield 

^ ^ 3(5.)ifle4°o'(f)g?(f).i . 6 X 10- . (46) 

2{GMr^)2 

for d = 5 km and rp = 25498 km. As can be seen, AP/^^^ » AP/^^-^~^'\ 
so that the impact of the direct tidal perturbations on the measurement of 
the gravitomagnetic clock effect on I can be neglected. 
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The cross-coupling between the static Earth oblateness and the Earth 
tides would induce the following perturbations on the mean longitude 

dl 21J2GMR'+^A,^ Y!p=o Et" oo FimpG,pq[{t - 2p) - mjkf^Hp 

— = n j-T cos 

dt 2r^+5f 

(47) 

Prom an inspection of the explicit expressions of the inclination functions 
(see Appendix A) and from the condition i — 2p = 0, which holds for even 
i, it turns out that the i = 2 part of the indirect tidal spectrum, which is 
the most powerful, has no effect on / for an equatorial orbit. More precisely, 
F201 does not vanish for i = 0°, but in this case, £ — 2p = m = 0. The other 
i = 2 inclination functions vanish for equatorial orbits. The same holds also 
for the i = 3 constituents for which i — 2p^0. It turns out that F^n ^ for 
equatorial orbits; in this case £ — 2p = m = 1. The other £ = 3 inclination 
functions vanish for i = 0° . 

3.3.3 The Impact of the Errors in the IncUnations 

Until now we have propagated the uncertainties in the orbit radius by assum- 
ing for the inclination and the eccentricity their nominal values i = e = 0. 
Now we wish to fix tq and propagate the errors in the inclination As 
pointed out before, the uncertainties in i affect the classical perturbative 
corrections to the Keplerian periods. The effects of the Earth oblateness 
are the most prominent; thus, we will investigate them in order to establish 
upper bounds in the admissible separation / between the two orbital planes. 

The radial component of the perturbing acceleration due to J2 for a 
generic orbit reads 

^ _ 3GMR(j^J2 [3cos2(a; + /) +3cos2i-3cos2icos2(a; + /) - l] . 

(48) 

Neglecting all terms of order 0{e^), ea. H48() in eq.© yields 

{G Mr 0)2 

Then, for the counter-orbiting satellites along identical (almost) circular 
orbits with inclinations i^"^^ = i and i^"-* = (180° — i) + I, with I /i -C 1, 
from cos^ « ~ 1 + i^, 

' (GMro)5 
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For ro = 25498 km and i = 0.01° , AP/-^'^ < AP^^ ' if / < 0.006° = 1 x lO"-^ 
rad. It is worth noting that the current technology does allow to obtain 
equatorial orbits tilted by 0.01° to the equator, e.g., for many geostationary 
satellites. Many orbital perturbations affect the inclination with shifts Ai 
which must be kept smaller than /. The static part of geopotential does 
not induce long-term perturbations on the inclinations. The tesseral and 
sectorial solid Earth and ocean tides induces long-period perturbations on i 
of the order of a few tens of mas. The general relativistic gravitoelectric de 
Sitter, or geodetic, precession [^5] induces a secular variation of the incli- 
nation of 84 mas yr~^. The inclination is sensitive to the non-gravitational 
perturbations, but the drag-free apparatus could resolve this problem. Suf- 
fices it to say that the uncancelled non-conservative perturbations on i for 
the LAGEOS satellites amount to a few tens of mas yr~^ [201 • It is also 
important to note that departures of i from their nominal values of the or- 
der of / ~ 0.0001° are feasible with the current technologies for, e.g., the 
planned GP-B mission [371 138j and the current GRACE mission^^ . It turns 
out that the requirements on the semimajor axis and the eccentricity are far 
less demanding than those on the inclination j34j . 

3.3.4 The N-Body Gravitational Perturbations 

Another source of perturbations on an Earth satellite's mean longitude is 
represented by the gravitational effect induced by the major bodies of the 
Solar System (Sun, Moon, Jupiter, other planets, the asteroids). Let us 
calculate the effects induced by some of them. 

The perturbative effect of the planet of mass m on the satellite of mass 
m is given by^^ [25j 



It turns out that the second term in ea. lfS^ does not induces secular per- 
turbations. 

^^http : / /www. csr .utexas . edu/grace/newsletter/2002/august 2002 .html 

^^Here we adopt the Lagrangian approach; TZ is the disturbing function and the (ap- 
proximate) equation for the rate of the mean longitude is 




(52) 



di_ _ 2 on 

dt na da 



(51) 
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After expressing the first term of eq. (|52|) in terms of the orbital elements 
of the satellite and the planet and averaging it over one period of the mean 
longitudes I and / one finds for the largest contribution which does not 
contain the terms of second order in the eccentricities and the inclinations 

(planets) ^ _47rGm ^ 

It is important to notice that eq. (|53|) does not depend on the sense of motion 
along the orbits, so that it cancels, in principle, in APi. Instead, the differ- 
ence d in the semimajor axes of the two satellites would induce an aliasing 
effect 

^^(planets) _ ISlTm dJ^ ^ ^^^^ 

The nominal values of eq . (|54|) for the Sun and the Moon amount to 1.178x 10~^ 
s and 2.565 x 10~^ s, respectively. However, the uncertainties in Gm yield 

5[APf'^°)]GMo = 7xl0-i^s, (55) 
5[AP/^°™^)]c™_„ = 6 X 10-11 s, (56) 

where we have used 5{GM)q = 8 x 10"^ m'^ s~^ [SO] and 5{Gm)Moon = 
1.2 X 10"^ m^ s-2 ISU. The errors induced by the uncertainties in d, a and 
GM^ turn out to be negligible. The corrections /^p^P^^'^^^^'^ Venus, Mars, 
Jupiter and Ceres, whose mass amounts to approximately one third of the 
total mass of asteroids TTiasteroids _ 2 3 x 10^'^ g, are negligible because their 
nominal values are < 10~^ s. This is particularly important for Jupiter, 
whose contribution would amount to A.P^'^^^^ = 1 x 10~^ s, because its 
sidereal orbital period, with respect to the Sun, amount to almost 11 years. 

It can be shown that the indirect effect induced on P*^^^ by the pertur- 
bations on a can be neglected. Indeed, there are no N-body secular pertur- 
bations on a and, most importantly^^, their amplitudes are proportional to 
factors like el'^^'e'l'^^l ginl'^^l i sinl'*'*! i', where /ii, /12, and /14 are integer 
numbers [HI] constrained by the conditions + I/i2| + l^sl + < 2 and, 
for the short-period perturbations, /ii + /i2 + ^3 + /14 / 0. 

3.4 The Impact of the Non-Gravitational Perturbations 

From eq.© it turns out that, in general, periodically time- varying radial 
accelerations with frequencies multiple of the orbital frequency do not affect 

"'^'^ After all, short-periodic terms in the planetary longitudes would result in semi-secular 
signatures on the time-scales of an Earth satellite. 
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the mean longitude over an orbital revolution, at least at order 0{e). 

The situation is different for radial accelerations which can be considered 
constant in time over an orbital revolution. Let us consider an acceleration 
of non-gravitational origin whose radial component Ar is constant over the 
timescale of one orbital revolution. From eq.Q it turns out that 

from which it follows 

AFr' = 14..(^)^«. (58) 

The maximum value of Aji which makes APj^^'^^ < APj^'^^^ is A^'^^ = 
6 X 10~^ cm s~'^ for d = 5 km and = 25498 km. It is not a too stringent 
constraint. For LAGEOS, which is completely passive, the largest non- 
gravitational acceleration, i.e. the direct solar radiation pressure, amounts 
to almost 4 x 10~^ cm s^^; a drag- free apparatus with not too stringent 
performances- well far from the 10~^^ cm s~'^ Hz~2 level of LISA-could 
meet in a relatively easy way such requirement. For 6d ~ 6rQ ^ 1 cm, 
6{GM) = 8 X 10^^ cm^ s"^ and Ar = 6 x 10"^ cm it turns out 

<5[AP/^^)], = Utt (^-^y Anx{6d) = 1x10-'^ s, (59) 
S[AP[''%, = 357rd(^)'ylRx(5ro) = 5xlO-i6g^ (gO) 

,5[APf^)]GM = 2l7Td(^-^y ARx5iGM) = 1x10-''' s. (61) 

These results show that the direct impact of the non-gravitational accel- 
erations on the measurement of the gravitomagnetic time shift on / could 
be made negligible with a drag-free apparatus of relatively modest perfor- 
mances, contrary to the indirect effects on the difference of the Keplerian 
periods which, as seen before, would require a much more effective drag-free 
cancellation in the frequency range [0, 10~^ Hz]. 

4 Conclusions 

In this paper we have examined the possibility of measuring the gravito- 
magnetic clock effect on the mean longitudes of a pair of counter-rotating 
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satellites following almost identical circular equatorial orbits in the gravi- 
tational field of Earth. While the gravitomagnetic signature depends only 
on the Earth parameters, the aliasing classical effects depend also on the 
orbital geometry of the satellites, so that it is possible to choose it suitably 
in order to reduce their impact on the measurement of the post-Newtonian 
effect. In this respect, our choice of a nominal value of 25498 km for the 
semimajor axis yields good results. 

From the point of view of the observational accuracy, one source of error 
would be the uncertainty with which the axes of the ITRS are known. How- 
ever, it turns out that the level of the gravitomagnetic effect can be reached 
over a sufficiently high number of orbital revolutions. 

In regard to the systematic errors, the major limiting factors are 

• The unavoidable imperfect cancellation of the Keplerian periods, which 
yields a 10^^ cm constraint in knowing the difference d between the 
semimajor ax;es of the satellites 

• The required accuracy in the knowledge of the inclinations i of the 
satellites in presence of not exactly equatorial orbits. For i = 0.01° 
the difference between the inclinations of the two orbital planes / = 

- should be less than 0.006°. 

A pair of drag-free (10^^ cm s^^ Hz^2 ) spacecraft endowed with an inter- 
satellite tracking apparatus might allow to meet the stringent requirements 
of such a mission whose realization seems to be very difficult although, per- 
haps, not completely impossible with the present-day or forthcoming space 
technologies. 
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